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ABSTRACT 

We show that the distribution of the sizes and temperatures of clusters can be used to con- 
strain cosmological models. The size-temperature (ST) distribution predicted in a flat Gaus- 
sian cluster-abundance-normalized ilo = 0.3 model agrees well with the fairly tight ST rela- 
tion observed. A larger power-spectrum amplitude erg would give rise to a larger scatter about 
the ST relation as would a larger value of f^o and/or long non-Gaussian high-density tails in 
the probability density function. For Gaussian initial conditions, the ST distribution suggests 
a constraint af^^^-^^ ~ 0.76. The ST relation is expected to get tighter at high redshifts. In 
the process, we derive a simple formula for the halo formation-redshift distribution for non- 
Gaussian models. We also suggest that the discrepancy between the naive zero-redshift ST 
relation and that observed may be due, at least in part, to the fact that lower-mass clusters 
form over a wider range of redshifts. An Appendix derives an equation for the formation- 
redshift distribution of halos. 
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1 INTRODUCTION 

Galaxy clusters are now being widely used as probes of cosmolog- 
ical and structure-formation models. For example, the abundance 
. of galaxy clusters has been used to constrain the amplitude erg of 
' the power spectrum and the nonrelativistic-matter density fio in 
models with an initially Gaussian distribution of density perturba- 
tions (Evrard 1989; Henry & Amaud 1991; Bahcall & Cen 1992; 
Bahcall & Cen 1993; Lilje 1992; Oukbir & Blanchard 1992; White, 
Efstathiou & Frenk 1993; Viana & Liddle 1996; Eke, Cole & Frenk 
1996; Viana & Liddle 1999), as well as in models with long non- 
Gaussian tails (Robinson, Gawiser & Silk 2000; hereafter RGS). 

In this paper, we show that the scatter in galaxy cluster scaling 
relations can be used to constrain cosmological and structure for- 
mation models. Specifically, we focus on the relatively small scatter 
of the relation between X-ray isophotal size and emission-weighted 
intracluster-medium mean temperature Tx demonstrated in Mohr 
& Evrard (1997; hereafter ME97). We illustrate how this scatter 
should depend on and f2o, and how it is affected by the intro- 
duction of a non-Gaussian distribution of perturbations with a long 
tail of high-density peaks. Our work on the size-temperature (ST) 
relation follows prior analytic work by Kitayama & Suto (1996) 
(although they focussed primarily on other cluster properties) and 
employs the framework for relating the ST relation to the underly- 
ing dark matter properties as discussed in Mohr et al. (2000; here- 
after MOO). 



The small scatter is heuristically expected if clusters form at 
rare high-density peaks in a Gaussian primordial distribution. Clus- 
ters that form earlier should be denser when they are first virialized 
and so they should have smaller radii for a given mass, or similarly, 
smaller radii for a given temperature. In this way, any dispersion 
in the formation redshifts for clusters of a given mass should yield 
a spread in the ST relation. If clusters come from rare Gaussian 
peaks, then the spread in formation redshifts should be small; given 
the rapidly dying Gaussian tails, it is unlikely that any cluster of a 
given mass observed today was formed at a redshift much earlier 
than the others. However, if the distribution had long non-Gaussian 
tails (as would be required to significantly boost the cluster abun- 
dance) or if clusters formed from peaks that were not quite so rare 
(e.g., > 2a rather than > Scr peaks), then clusters of a given mass 
observed today should have had a much broader distribution of for- 
mation redshifts (see Fig. [l]) and thus a much broader distribution 
of sizes (for a given mass or temperature). 

We quantify these arguments using a spherical-top-hat- 
collapse model to relate the virial radius and temperature of a clus- 
ter to its mass and formation redshift. We use the formation-redshift 
distribution for Gaussian perturbations from Sasaki (1994), and we 
generalize it for an arbitrary initial density distribution (the deriva- 
tion is presented in an Appendix). We use a Monte Carlo approach 
to simulate the ST relation for a variety of parameters, and illustrate 
in particular how it depends on ag, f2o, and G, the non-Gaussian 
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Figure 1. The solid curve shows a Gaussian distribution P(y) with unit 
variance, while the broken curve shows a non-Gaussian distribution with 
the same variance but 10 times as many peaks with y > 3. This illustrates 
(a) how the cluster abundance can be dramatically enhanced with long non- 
Gaussian tails (since clusters form from rare peaks); and (b) that the disper- 
sion of y for y > 3 is much larger for the non-Gaussian distribution than it 
is for the Gaussian distribution, and this will lead to a larger scatter in the 
formation redshifts and sizes of clusters of a given mass. 

multiplicative excess of > 3cr peaks introduced by RGS. Our main 
results are (a) the predicted scatter in the ST relation for Gaussian 
initial conditions and favored cosmological parameters is found to 
be fairly consistent with that observed; (b) G > 5 greatly over- 
predicts the scatter; (c) the scatter for the non-Gaussian initial con- 
ditions required to make the cluster abundance consistent with an 
Einstein-de-Sitter Universe (EdS) is also much larger than that ob- 
served. Joint constraints from the cluster abundance and the ST re- 
lation on (jg, G, and Q,o are discussed. We show how the ST re- 
lation should be altered for clusters at intermediate and high red- 
shifts. In the process we show that, because lower-mass clusters 
form over a larger range of redshifts than higher-mass clusters, the 
expected ST relation is steeper (and therefore more consistent with 
the observed relation) than the naive expectation detailed in MOO. 
In the final Section, we make some brief connections to the X-ray 
mass-temperature relation and to the redshift evolution of the clus- 
ter abundance. 



2 INGREDIENTS 

2.1 Spherical-Collapse Model 

We use the relations of Kitayama & Suto (1996) to relate the cluster 
virial radius and virial temperature at formation time, i?vir and T, 
to the mass M and formation redshift Zf (defined to be the redshift 
at which the cluster collapses). Fig. ^ shows how this model assigns 
masses and formation redshifts to clusters of given temperatures 
and sizes assuming that J?vir oc Rs . 

It is possible to connect more rigorously these cluster dark- 
matter properties with the observable intracluster medium (ICM) 
properties in a manner similar to that outlined in MOO. Specifically, 
we assume that Tx is the virial temperature (e.g. Evrard, Metzler 
& Navarro 1996; Frenk et al. 2000; Bower et al. 2000). We trans- 
form from the virial radius at z/ to the X-ray isophotal size R us- 
ing R oc R^ri^ f'icM' where ficM is the ICM mass fraction (MOO 
eqns. (8) and (10)). The dependence on ficM should be included 
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Figure 2. Mass and formation-redshift contours in the size-temperature 
plane for f!o = 0.3 and h = 0.65 obtained from the spherical-top-hat 
model of gravitational collapse discussed in the text. It is clear from the fig- 
ure that a narrow (broad) spread in the formation redshift will yield a tight 
(broad) ST relation. For larger Oq > the Zj = Q contour remains the same, 
but the spacing between ejofi\-Zf contours increases. 

because variations in ficM with mass are observed (e.g. Mohr, 
Mathiesen & Evrard 1999; David, Jones & Forman 1995) and could 
alter the slope of the ST relation. In the following analysis, we as- 
sume ficM OC TjJ;'^'*; however the results of Figs. 4 and 5 are no 
more than very weakly dependent on the fiCM functional form. 

We normalize the simulated ST relation to the observations by 
fixing the constant of proportionality so that no observed cluster in 
the local sample lies above the = line (see Figs 2 and 3). 

2.2 Distribution of Halo Masses 

Numerical simulations tell us that the Press-Schechter (PS) ap- 
proach (Press & Schechter 1974) provides a reasonable approxi- 
mation for the abundance of cluster size halos of a given mass at 
any given epoch for Gaussian initial conditions (e.g. Lacey & Cole 
1994, Gross et al. 1998, Lee & Shandarin 1999), and for a few non- 
Gaussian initial conditions that have been explored with simula- 
tions (Robinson & Baker 1999). In the PS approach the number per 
comoving volume of halos with masses between M and M + dM 
at redshift z is (e.g., Lucchin & Matarrese 1988; RGS), 

dM M V 'j dM 

where pt is the background density, P{y) is the primordial proba- 
bility distribution function normalized to unit variance. The argu- 
ment y = 5{z)/aM, and 5{z) — 5{z)/D{z) where 5c{z) is the 
critical overdensity for collapse (see Kitayama & Suto 1996 for ac- 
curate analytic fits), and D{z) is the linear-theory growth factor. 
Here, gm is the current root- variance of spheres that enclose an 
average mass M, and f ~ P{y) dy. 

2.3 Distribution of Formation Redshifts 

The objects of mass M observed at some given redshift Zo under- 
went collapse at a variety of formation redshifts zj > Zo. Sasaki 
(1994) has shown how the PS formalism leads to an expression for 
the formation-redshift distribution under the assumption of Gaus- 
sian initial conditions and that the merger rate has no characteristic 
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Figure 3. (a) ST distribution for LCDM and erg = 0.99 and Gaussian initial conditions. Eacli dot represents a simulated cluster, while the diamonds are data 
from MOO. The line shows the ST relation expected for clusters that form today, at redshift 2 = 0. (b) shows the same except that here we use erg = 1.5. (c) 
shows the same as in (a) but with the non-Gaussian distribution of RGS with G = 10. 



mass scale. His derivation can be generalized in a straightforward 
fashion to arbitrary P{y) . Doing so (see Appendix), we find the dis- 
tribution (normalized to unity) of formation redshifts z f for halos 
of mass M observed at redshift Zo to be. 



'(y(Af, 



(2) 



where P'(y) = dP/dy. Lacey & Cole (1993,1994) have presented 
an alternative, but somewhat more complicated, formation-redshift 
distribution that improves upon Sasaki's assumption of self-similar 
merging. We will leave the implementation of this alternative dis- 
tribution and a discussion of the formalism introduced by Percival 
et al. (2000), to future work, but note that our preliminary investiga- 
tions, as well as previous results (Viana & Liddle 1996; Buchalter 
2000), indicate that the predictions of these models do not differ 
considerably for cluster-mass halos. 



2.4 Preliminary Estimates 

It is straightforward to roughly estimate the effects of non-Gaussian 
tails on the ST-relation scatter. For a rapidly dying distribution 
P{y), the controlling factor in dn/dzf will be P'{y)- For a Gaus- 
sian P(y), the root-variance of y is 0.282 for the distribution P'(y) 
for values of y > 3, and the mean value of y is 3.30. For an EdS 
model, y — 1.69(1 + Zf)/aM, and (1 + z)~^ oc i?vir. Thus, 
an/R = {4:/3){aR^,JR,ir) ^ i4:/3)(ay/y) = 0.113 for a Gaus- 
sian distribution, in surprisingly good agreement with the estimate 
of the intrinsic scatter of 10% in the ST relation (ME97). For the 
RGS distribution with G — 10, the root-variance is 0.896 and the 
mean value of y is 3.87 leading to CTR^^^/Rvir — 0.31, more than 
twice the observed scatter. Below we will quantify this far more 
precisely. 



3 RESULTS 

For any given fio, ""8, and G, we perform a Monte Carlo realization 
of 400 clusters with the mass and formation-redshift distributions 
given above. We then assign to each of these clusters a size and 
temperature as outlined in Section '. 



2.1 



The ME97 sample to which 
we compare our calculations is a flux-limited sample. Within this 



sample, the probability of finding a cluster of luminosity Lx goes 
as Lx^, and Lx is observed to go as roughly T^'^ to (David et 
al. 1993; Amaud & Evrard 1999), so the flux limit is essentially a 
virial-temperature weighting of T^-'^^ to T*'^. We thus subject our 
simulated population of clusters to a T'^-'^^ weighting; our results 
are not significantly altered for the steeper weighting T*'^. 

Fig. ^a) shows the results of our Monte Carlo for a flat 
flo = 0.3 model (LCDM) with the value ag = 0.99 inferred 
from the cluster abundance (Viana & Liddle 1999) and a Gaus- 
sian distribution. The data points from ME97 are overlaid. We used 
a Hubble parameter h = 0.65, but the results are essentially un- 
altered for different plausible values of h. Fig. ^(b) illustrates that 
the scatter in the ST relation is increased if the power-spectrum 
normalization is higher. In this case, clusters are not quite as rare, 
and they form over a larger range of redshifts. Fig. ^(c) shows how 
the scatter is increased as the abundance of high-density peaks is in- 
creased. In this case, clusters observed today are also formed over a 
broader range of redshifts. At this point, we note the apparent sim- 
ilarity between the predictions of the ST distribution of the cluster- 
abundance-normalized Gaussian LCDM model and the data; the 
scatter about the ST relation would be broadened considerably with 
a higher as or with a highly non-Gaussian model. 

To make these arguments more quantitative as well as survey 
a larger range of parameters, we have simulated ST relations for a 
variety of models in the crg-G parameter space for both EdS and 
LCDM models and then used a 2D Kolmogorov-Smirnov (KS) test 
(Peacock 1983; Press et al. 1996) to compare these with the data. 
Fig. ^ shows the resulting contours of constant KS significance 
levels for both f2o = 0.3 and fio = 1. The results suggest that 
the Gaussian cluster-abundance-normalized (ug — 0.99) LCDM 
model provides a good fit to the data. 

We heuristically expect that the dependence of the ST scatter 
on cosmological parameters/models should be similar to that of the 
cluster abundance; if the peaks that give rise to clusters are rare, 
we expect little scatter and vice versa if clusters are more com- 
mon. The contours of fixed cluster abundance in Fig. ^ indicate 
that this is qualitatively correct. We obtain these curves by using a 
cluster abundance n(> 6.2keV,2 = 0.05) = 1.53 x IQ^O '^^ x 
10~^ Mpc~^ (Viana & Liddle 1999) and integrating eqn. (§) up 
from the mass associated with a temperature 6.3 keV and a forma- 
tion redshift 2/ = 0. However, the detailed results also seem to 
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Figure 4. The heavy sohd curves show confidence levels suggested by 
the ST data in the o-g-C parameter space for the LCDM model, and the 
light curves show the same for the EdS model. The dashed curve shows 
the contour suggested by the central value [n(> 6.2keV,2: = 0.05) = 
1.53 X 10±0-1S X 10"'^ h^} of the local cluster abundance for 

f!o = 0.3, while the dotted curves indicate contours for the upper and 
lower observational limits to the cluster abundance. 



indicate that if Q.q is fixed, the ST distribution and cluster abun- 
dance can be used in tandem to break the degeneracy between G 
and (T8. In fact, combining the two constraints already seems to rule 
out large deviations from Gaussianity. 

Fig. ^ shows the regions of as-^o parameter space preferred 
by the ST relation, as well as the curve in this parameter space 
suggested by the cluster abundance. For fixed ag, the ST scat- 
ter increases as f2o increases. At first, this might seem discrepant 
with the well-known result that the range of formation redshifts 
is narrower for larger fio for cluster-abundance-normalized mod- 
els. However, this narrowing of the formation-redshift distribution 
with increasing S7o is not quite as dramatic if we fix erg instead 
of the cluster abundance. More importantly, the spherical-top-hat- 
collapse dynamics leads to a broader spacing between the equi-z/ 
contours in Fig. ^, and this is responsible for increasing the ST scat- 
ter as 0,0 is increased with fixed ag; in other words, the relationship 
between R and T evolves more rapidly with redshift in higher Q.o 
models. 

From the results in Fig. ^, we can approximate an ST con- 
straint, (Tg = 0.76 Q.q'^''^^ , as compared with the cluster-abundance 
constraint, as = 0.56 ° (Viana & Liddle 1999). The region 
of overlap between the cluster-abundance constraint and the ST re- 
lation lies at low values of f2o, low values of non-Gaussianity, and 
slightly higher values of ag . 

3.1 An Einstein-de-Sitter Universe? 

RGS were able to identify for an EdS model, a region in the ag-G 
parameter space near ag = 0.4 and G = 10 in which the predicted 
cluster abundance was found to agree with that observed. Fig. ^ 
shows that these parameter choices predict far too much scatter in 
the ST relation. Allowing for additional sources of scatter in this 
simulated ST relation would only increase the discrepancy between 
the model and the observations. 



Figure 5. The heavy solid curves show likelihood contours suggested by the 
ST data for Gaussian initial conditions in the crg-fio parameter space. The 
dot-dashed curve shows the contour preferred by the local cluster abundance 
as suggested by Viana & Liddle (1999), while the dashed curve shows the 
fit to our ST constraint. 
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Figure 6. The ST distribution for Qo = 1 with erg = 0.4 and G = 10, 
one of the combinations of parameter values that yield the correct cluster 
abundance for an EdS Universe. The predicted scatter in the ST relation is 
considerably larger than that observed. 

3.2 High and Intermediate Redshift Results 

Clusters that exist at higher redshifts must form from even higher- 
density peaks than those today. Thus, in a Gaussian model, the scat- 
ter in their formation redshifts and thus in their sizes should be even 
smaller. This is illustrated in Fig. ^. The canonical-model predic- 
tions shown in Fig. Q(a) for z ~ 0.3 seem to be in relatively good 
agreement with the cluster sample observed so far. Fig. ^(b) shows 
that the scatter in the ST relation for the canonical model should be 
very small. Even though the sample of such high-redshift clusters 
is expected to be small, the predicted scatter is so small that mea- 
surement of the sizes of only a handful of clusters could put strong 
constraints on different sources of scatter (e.g., non-Gaussianity, 
measurement uncertainties, mergers, galaxy feedback, etc.). 

3.3 A size-temperature anomaly? 

The ST relation of the low-redshift X-ray flux limited cluster sam- 
ple has a slope of m ~ 1, which is considerably steeper than the 
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Figure 7. The ST relation for the LCDM model with as = 0.99 and 
Gaussian perturbations for (a) clusters at 2 = 0.3 and (b) 2 = 1. This 
illustrates how the scatter should decrease with redshift. The data from 
Mohi' et al. (2000) (that have median redshift 0.3, but cover the range 
0.19 < 2 < 0.54) are shown in (a). 



m = 2/3 slope expected in a model where all clusters are assumed 
to have formed at the redshift of observation (ME97, MOO). ME97 
suggest that a possible explanation for this steeper than expected 
scaling relation is provided by galaxy feedback. Fig. ^ illustrates 
that the discrepancy may be due only, or at least in part, to the 
fact that lower-mass clusters form over a broader range of redshifts, 
and thus will in general have smaller sizes than they would if they 
all formed very recently. Visual inspection of Fig. suggests that 
this is a plausible explanation, especially when the small-number 
statistics of the observational sample are taken into account. More- 
over, the relatively strong dependence of the ST scatter on erg indi- 
cate that better agreement than shown in Fig. |3| could be obtained 
with a slightly different value of o-g and/or Q,q (cf.. Fig. The 
apparent disagreement with ME97's feedback-free numerical sim- 
ulations, which show an ST scaling close to the naive scaling (but 
still steeper in 3 of the 4 cosmologies tested), may have been due 
to this (78 and fio dependence and/or the relatively small-number 
statistics of their simulations sample. Thus, the apparent deviation 
of the ST relation slope from the m = 2/3 expectation is not quite 
so anomalous. 



4 DISCUSSION 

We have calculated the ST distribution of clusters with a simple an- 
alytic model and focussed in particular on the dependence on the 
power- spectrum amplitude (erg) and the degree of non-Gaussianity 
(G). We find a fairly sensitive dependence of the ST relation scat- 
ter on these two parameters. Thus, the tightness of the ST rela- 
tion can be used to place valuable constraints on these parame- 
ters, as well as on other cosmological parameters. The canonical 
cluster-abundance-normalized Q.q = 0.3 model predicts an ST re- 
lation consistent with that observed, but a erg much larger or smaller 
would be inconsistent as would a non-Gaussian model that predicts 



a significant excess of > 3(t high-density peaks.^ The constraints 
to erg, G, and fio that arise from the ST distribution should be qual- 
itatively similar to those from cluster abundances, but our prelim- 
inary calculations suggest that they may be sufficiently different 
to provide complementary constraints. The ST relation should be- 
come increasingly tight at larger redshifts. Our results also suggest 
that the discrepancy between the naive 2 = ST relation and the 
observed ST relation may be due, at least in part, to the fact that 
lower-mass clusters observed today have formed over a larger range 
of redshifts than higher mass clusters. 

The fact that lower-mass clusters tend to form over a broader 
redshift range than higher-mass clusters will also tend to steepen 
the Mvir-T relation beyond the self-similar expectation of m = 
2/3. Numerical simulations of structure formation within models 
with non-Gaussian initial conditions or \ow-Qo open models ought 
to exhibit this effect. The OCDM256 portion of Fig 4. in Bryan & 
Norman (1998) indicates that low-mass clusters fall systematically 
below the best fit AI^h-T relation, consistent with our expectation. 
It should be emphasized that, in this particular study, the low-mass 
systems are composed of far fewer particles than the high-mass sys- 
tems, providing another plausible explanation for structural differ- 
ences. Further work to investigate departures from self similarity 
in the cluster population which naturally arise from the spread in 
formation epochs is clearly required. 

Since the overdensity-peak amplitude at which a cluster can 
form increases at higher redshift, the redshift evolution of the clus- 
ter abundance depends on the shape of the primordial density distri- 
bution function at high peaks just as the ST scatter does. Thus, if fio 
is fixed, it should be possible to reconstruct the cluster-abundance 
evolution from the scatter in the ST relation for local clusters. 

Although we have used cluster sizes inferred from X-rays to 
compare with theoretical calculations, the same could be done for 
the sizes of clusters measured via the Sunyaev-Zeldovich effect, 
either with or without redshift information (Kamionkowski et al. 
2000). Of course, there will invariably be some cluster-formation 
physics that our current analysis has left out, and numerical sim- 
ulations may have an advantage in this regard. Note that the only 
source of scatter in our simulated ST relations is the range of forma- 
tion epoch, whereas other sources of stochasticity (e.g., measure- 
ment uncertainties, mergers, galaxy feedback, etc.) might increase 
the scatter. However, with our analytic approach, we are able to 
rule out models that overpredict the scatter; we can sift far more 
rapidly through a variety of cosmological models and parameters, 
study the dependence of the ST distributions on these models and 
parameters, and gain some intuitive feel for how the results arise. 
By doing so, we hope to have established that cluster sizes can pro- 
vide a valuable new probe of cosmological models. 
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APPENDIX 

Taking the derivative of eqn. ^ with respect to redshift, we ob- 
tain (hereafter we do not explicitly show the {M, z) dependencies 
where they are obvious) 



dMAz 



fPb 
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dP_dydy_ d dy 

dy dz dM dz dM 
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1 d& do 
~dzdM 



dy a 
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(4) 



(5) 



AMAz AMAz 
In the last line we have equated the total rate of change to the dif- 
ference between a formation rate and a destruction rate (the latter 
being due to objects merging to form larger objects). These can be 
expressed as 

nn. 

(6) 



d nfo 



AMAz 



^Q{M,M';z)dM', 



where Q{M, M'; z) is the probability that an object of mass M' 
is one of the merging components when an object of mass M 
forms and Mmin is introduced to prevent the integral from di- 
verging, and d^ridcst/dMdz = <j){M, z){An/dM), where the 
function </)(M, z) can be interpreted as the destruction rate per 
bound object. Sasaki assumes that <j>{M, z) can be expressed as 
(j){M,z) — M"<j){z) (implying that the destruction rate has no 
characteristic mass scale). Using eqn. (^), we can write 

- -d^Ti/dAfdjz + d^nform/dMdz 



An/AMM° 



(7) 



Since the left-hand side of this equation depends only upon z the 
right-hand side must be independent of M and so may be evaluated 
at a very small mass Afmin. Since the formation rate is zero at Mmm 
(see eqn. H), this leaves 



7. X _ d^n/dAMz(Mmin,z) 
An/AM{M^i^,z)M^,J 

Substituting eqns. (|l|) and (^ into this expression gives 

1 9P[y(Afn,in)] 5 



(8) 
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dy 



+ 1 



(9) 



For a hierarchical clustering model, limAf-.o (t^(A/) — oo, so if 
we take the limit Afmin — > eqn. (^ will be or oo unless q = 0, 
forcing the choice a — upon us such that (f>(z) — (1/S){AS /Az). 
Substituting this expression and eqn. (^ into eqn. (Q), we find that 
the formation rate is given by 



d^riform _ fpb 1 dS da dP S 
AMAz ~ ~W^dzdM'dya' 



(10) 



This is the rate of formation of bound objects of mass M and red- 
shift z, but we wish to know what fraction of these objects will 
survive until the redshift of observation. Using our definition of 
(piz) the number of objects of mass M which formed at Zf must 



evolve with redshift as AN/Az — <j>{M, z)N such that the frac- 
tion remaining by Zo{< Zf) is f{zi,Zo) = exp (j>{z)Az = 
S{zo)/S(zi). The number of objects of mass A/, which formed at 
redshift Z{ and which survive until redshift Zo is given by the prod- 
uct of this expression and eqn. (p^; i.e.. 



d^n _ fpb S{zo) dy. , dy . ,dP. 
AMAz M 5{zf) dz ^ " dM ^^'^ dy 



(11) 



Noting that {dy /dM),J{dy /dM),^ = Mzi) / 5{zo), we obtain 
our final result, eqn. (Q), by dividing eqn. ( plj ) by eqn. (|l|). 
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